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Abstract 

Background: Myricaria consists of about twelve-thirteen species and occurs in Eurasian North Temperate zone, most species 
in the Qinghai-Tibet Plateau (QTP) and adjacent areas. 

Methodology/Principal Findings: Twelve species of Myricaria plus two other genera Tamarix and Reaumuria in 
Tamaricaceae, were sampled, and four markers, ITS, rpsl6, psbB-psbH, and trnL-trnf were sequenced. The relaxed Bayesian 
molecular clock BEAST method was used to perform phylogenetic analysis and molecular dating, and Diva, S-Diva, and 
maximum likelihood Lagrange were used to estimate the ancestral area. The results indicated that Myricaria could be 
divided into four phylogenetic clades, which correspond to four sections within the genus, of them two are newly described 
in this paper. The crown age of Myricaria was dated to early Miocene ca. 20 Ma, at the probable early uplifting time of the 
Himalayas. The Himalayas were also shown as the center of origin for Myricaria from the optimization of ancestral 
distribution. Migration and dispersal of Myricaria were indicated to have taken place along the Asian Mountains, including 
the Himalayas, Kunlun, Altun, Hendukosh, Tianshan, Altai, and Caucasus etc., westward to Europe, eastward to Central 
China, and northward to the Mongolian Plateau. 

Conclusions/Significance: Myricaria spatiotemporal evolution presented here, especially the Himalayan origin at early 
Miocene ca. 20 Ma, and then migrated westward and eastward along the Asian mountains, offers a significant evolutionary 
case for QTP and Central Asian biogeography. 



Citation: Zhang M-L, Meng H-H, Zhang H-X, Vyacheslav BV, Sanderson SC (2014) Himalayan Origin and Evolution of Myricaria (Tamaricaeae) in the Neogene. PLoS 
ONE 9(6): e97582. doi:10.1371/journal.pone.0097582 

Editor: Michael Hofreiter, University of York, United Kingdom 

Received April 19, 2013; Accepted April 22, 2014; Published June 6, 2014 

Copyright: © 2014 Zhang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This study was funded by China National Key Basic Research Program (No. 2012FY1 1 1500, 2014CB954201), CAS Important Direction for Knowledge 
Innovation Project (No. KZCX2-EW-305), and the Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences. The funders had no role in study 
design, data collection and analysis, decision to publish, or preparation of the manuscript. 

Competing interests: The authors have declared that no competing interests exist. 

* E-mail: zhangml@ibcas.ac.cn 



introduction 

The Tamaricaceae contains about eighty species [1] and four 
genera: Tamarix, Myricaria, Reumuria, and Hololachna [2]. This 
family, and Frankeniaceae, are defined as the salt-gland anatom- 
ical lineage [3] . Myricaria consists of about twelve - thirteen species 
[4-7] and occurs in Northern Temperate zone of Eurasia, mainly 
along the Asian mountains. There are eight species in Himalayas, 
many are endemic, thus forming a center of diversity for Myricaria 
(see Figure 1). 

Desvaux (1825) established the genus Myricaria and Niedenzu 
[8] presented the first classification. Zhang & Zhang [4] studied 
Myricaria in China and recognized ten species; they presumed the 
Himalayas to be the center of origin based mainly on the 
distribution of species. Gorschkova [7] described six species 
belonging to two sections in the flora of the former USSR. 

Another issue relevant to Myricaria systematics is the species 
Myricaria elegans. Ovezinrlikov & Kinzikaeva [9] erected the genus 
Myrtama based on this species but it caused some controversy. 
Zhang et al. [10] used ITS sequence data to study the relationships 
within Tamaricaceae and regarded Myrtama as an intermediate 



genus between Myricaria and Tamarix [11]. After sampling four 
species from Myricaria and sequencing ITS, rhcl^, and tRNAs Ser 
(GCU) and Gly (UCC), Gaskin et al. [1] found Myrtama and 
Hololachna to be distinct within Tamaricaceae, as did Zhang et al. 
[10]. However, based on additional sequence data, Hua et al. [12] 
and Wang et al. [1 3] confrrmed that Myrtama should be included in 
Myricaria. Sampling ten species of Myricaria and sequencing 
cpDNA p.^hA-trnli and the rpLl6 intron, Liu et al. [14—15] 
investigated the species-level phylogeographical patterns of Myr- 
icaria in western China as well as the origin oiM. laxiflora, a unique 
subtropical species of conservation concern from the Three Gorges 
of the Yangtze River in Sichuan and Hubei provinces. The 
Himalayas were proposed as the center of origin of Myricaria by 
Liu et al. [14] with the estimated age of origin 1.46-2.30 Ma. 

Closely associated with the distribution pattern oi Myricaria and 
related taxa, the QTP and Himalayan uplift during the Neogene 
are hypothesized to be a major influence on organism evolution in 
Asia (e.g. [16-17]). FoUowing collision of the Indian and Eurasian 
continents at ca. 50 Ma, the altitude and range of the QTP near 
the Oligocene-Miocene boundary became sufficient to trigger a 
reorganization of the Asian climate, as evidenced by the beginning 
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Figure 1. Distribution of Myricaria species (a,b), tKie information obtained from floras and herbaria, mainly in China. Geographical 
division of QTP eastern and western portions is shown in a red broken line (c). The Himalayan origin and dispersal routes along the Asian mountains 
are illustrated in arrows (a). 
doi:1 0.1 371 /journal.pone.0097582.g001 



of loess deposition in the Chinese Loess Plateau and the Junggar 
Basin [18-20]. Some evidence confirms that the central areas of 
the QTP were raised to present altitudes by that time [16,21-22] 
and uplift of the Himalayas may have also begun at that time [22- 
23]. Uplift of peripheral portions of the plateau has continued at 
various intervals [24—28]. A major uplift of QTP is often suggested 
to have occurred at 8 Ma, which also coupled with global cooling, 
even though Molnar [29] considered this uplift evidence to be 
inconclusive. Uplift of the QTP and global climate cooling and 
aridification [27] have been suggested causes for the evolution of 
many organisms [30-36]. As these studies have shown, rapid 
diversification of lineages in the QTP resulted in the migration of 
some species into other temperate regions, such as Central Asia, 
the Arctic, the Mediterranean (Caucasus-Alps) and southern Asia. 
Of these, connections between the QTP and adjacent arid, more 
northern areas can often be discerned, for example in recent 
studies on Hippophae rhamnoides (Elaeagnaceae) [37], Caragana 
(Fabaceae) [38], and Astragalus (Fabaceae) [39]. Linkage of the 
QTP and Africa and/ or the Mediterranean is illustrated by Begonia 
[40] and Uvaria (Annonaceae) [41]. An example linking the QTP 
and Southeast Asia is Paini (Anura: Dicroglossidae) [42] . 

The origin oi Myricaria has been associated with the QTP and 
Himalayas but justification has been weak. Zhang & Zhang [4] 
presumed that Myricaria originated from the Himalayas, only 
based on species distribution of the genus, whereas same opinion 
conducted by Liu et al. [14] from a phylogeography. Here we 
attempt to examine the origin and evolution of this genus and link 
it to the Himalayan uplift to explain the causes of its evolutionary 
patterns. In addition, the classification and distribution oi Myricaria 
are examined using molecular phylogeny and biogeography. 

Materials and Methods 

Taxa sampled 

Twelve species (seventeen samples) of Myricaria plus seven 
species from the outgroups Tamarix and Reummuria served as 
sources of DNA material (Table 1). The herbaria utilized in China 
were as follows: HNWT (Northwest Institute of Plateau Biology, 
Chinese Academy of Sciences (CAS), Xining, Qinghai); SHI 
(Shihezhi University, Shihezhi, Xinjiang); and XJBI (Xinjiang 
Institute of Ecology and Geography, CAS, Urumqi, Xinjiang), as 
well as the LE (Komarov Botanical Institute, Russian Academy of 
Sciences, St. Petersburg, Russia). 



DNA sequencing 

Total genomic DNA was extracted using the CTAB method 
[43]. The polymerase chain reaction (PCR) was used for 
amplification of double stranded DNA. The 25 |J,1 reaction system 
contained 0.25 |a.l of Ex Tag, 2.5 |J,1 of lOx Ex Taq buffer (Mg^"^ 
concentration of 25 mM), 2.0 |a.l of dNTP mix (2.5 mM concen- 
tration for each dNTP), 1 |j,l of the forward and reverse primers at 
5 umol/|J.l, and 0.5 |il of template DNA. The following primers 
were used: for ITS: ITSl-/(5'-AGA ACT COT AAC AAG GTT 
TCC GTA GC-3') and ITS4-r (5'-TCC TCC GCT TAT TGA 
TAT GC-3'), for tmL-trnF: tniL-f {5' -GGA AAT CGG TAG ACG 
CTA CG-3') and trnF-r (5'-ATT TGA ACT GGT GAC ACG AG 
-3'), for the intron of rpsl6: »^jl6-/ (5'-GTG GTA GAA AGC 
AAC GTG CGA CTT-3'), and for rpsl6-r (5'-TCG GGA TCG 
AAC ATC AAT TGC AAC-3') [44]; and the intergenic spacer 
psbB-psbH: psbB-psbH-r (5'-TTCAACAGTTTGTGTAGCCA-3') 
and /)jAB-/;jiH-/(5'-AGATGTTTTTGCTGGTATTGA-3') [45]. 

The protocol for amplification consisted of an initial hotstart at 
95°C for 2 min, followed by 30 cycles of denaturation at 94°C for 
30 s, annealing at 52°C for 30 s, extension at 72°C for 90 s, and a 
final extension at 72°C for 10 min. PCR products were purified 
using the PEG precipitation procedure [46]. These were 
sequenced using an ABI Prism 3770 Genetic Analyzer (Shanghai 
Sanggon Biological Engineering Technology & Service, Shanghai, 
China). Sequences were aligned using CLUSTAL X software [47], 
and then adjusted by hand. All gaps were treated as missing 
characters. Finally, a combined dataset consisting of four 
sequences of ITS and three cpDNA frnL-taF, rpsl6, and psbB- 
psbH, was prepared for phylogenetic analysis. 

Phylogenetic analysis and divergence time estimate 

The sequence dataset from twelve species (seventeen samples) of 
Myricaria plus seven species of Tamarix and Reummuria yielded 3202 
aligned nucleotide characters from four genes: ITS, taL-taF, 
rpsl&, and psbB-psb\i. The incongruence length difference (ILD) 
test of the four gene datasets was carried out using PAUP* [48], to 
assess potential conflicts between the different DNA fragments. 
This test was implemented with 100 partition-homogeneity test 
replicates, using a heuristic search option with simple addition of 
taxa, TBR branch swapping and MaxTrees set to 1000. 0.222 of 
incongruence length difference (ILD) tests [48] showed that the 
four gene datasets were not incongruent. 
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Table 1. List of Sampled Taxa, Vouchers and Genebank Accession Numbers. 





Species 


Voucher 


Source 


ITS 


trnL-trnF 


rpS^6 


psbB-psbH 


Myricaria Desv. 


M. alopecuroides Schrenk. 


P. Yan 3650 (SHI) 


Tashikurgan, Xinjiang, China, alt. 3650m 


KJ729654 


KJ729806 


KJ729756 


KJ729705 


M. alopecuroides Schrenk. 1 


Tibet-Xinjiang Exp. Team 1034 
(HNWP) 


Sukeplya, Yecheng, Xinjiang, China, alt. 
2800m 


KJ 808603 


KJ808634 


KJ 8086 19 


KJ808649 


M. bracteata Royle 


Y.H. Wu 36461 (HNWP) 


Nuomuhong, Dulan, Qinghai, China, alt. 
2840m 


KJ729655 


KJ729807 


KJ729757 


KJ729706 


M. elegans Royle 1 


P. Yan 3999 (SHI) 


Bandir, Tashekurgan, Xinjiang, China, alt. 
3000m 


KJ 808604 


KJ808635 


KJ 808620 


KJ808650 


M. elegans Royle 2 


P. Yan 7178 (SHI) 


Mazhaxi, Yecheng, Xinjiang, China, alt. 
3600m 


KJ 808605 


KJ808636 


KJ 808621 


KJ808651 


M. elegans Royle 3 


P. Yan 7378 (SHI) 


Ritu, Tibet, China, alt. 4600m 


KJ 808606 


KJ808637 


KJ808622 


KJ808652 


M. germenica (L.) Desv. 


1.0. Baitulin, Aralbaiev s.n. (LE) 


Zajsanskaya depression, E. Kazakistan 


KJ 808607 


KJ808638 






M. laxiflora (Franch.) P.Y. Zhang 
et Y.J. Zhang 1 


Wuhan Bot Card 


Wuhan Bot Card, Hubei, China 


KJ 808608 


KJ808639 


KJ808623 


KJ808653 


M. laxiflora (Franch.) P.Y. Zhang 
et Y.J. Zhang 2 


Wuhan Bot Card 


Wuhan Bot Card, Hubei, China 


KJ 808609 


KJ808640 


KJ 808624 


KJ808654 


M. paniculata P.Y. Zhang 
et Y.J. Zhang 


B.Z. Guo; W.Y. Wang 21930 
(HNWP) 


Linzhi, Tibet, China, alt. 2000m 


KJ 8086 10 






KJ808655 


M. platyphylla Maxim. 


Z.Y. Yang; LM. Ke 5711 (XJBI) 


Houxia, Urumqi, Xinjiang, China 


KJ80861 1 


KJ808641 


KJ808625 


KJ808656 


M. prostrata Hook.f. et Thomson 
ex Benth. et Hook.f. 


P. Yan 7242 (SHI) 


Hechakou, Hetian, Xinjiang, China, alt. 
5000m 


KJ 8086 12 


KJ808642 


KJ808626 


KJ808657 


M. pulcherrima Batalin 


LM. Ke 121 (XJBI) 


Ermuchang, Shaya, Xinjiang, China alt. 
4350m 


KJ 8086 13 


KJ808643 


KJ808627 


KJ808658 


M. rosea W.W. Sm. 


R.F. Huang G89-485 (HNWP) 


Milinpaiqu, Tibet, China, alt. 4530m 


KJ 8086 14 




KJ808628 


KJ808659 


M. squamosa Desv. 


P. Yan 4002 (SHI) 


Bandir, Tashekurgan, Xinjiang, China, alt. 
3002 m 


KJ729658 


KJ729810 


KJ729760 


KJ729709 


M. squamosa Desv. 1 


Y.H. Wu 3077 (HNWP) 


Beishan, Huzhu, Qinghai, China, alt. 
2700m 


KJ 8086 15 


KJ808644 


KJ808629 


KJ808660 


M. wardii C.Marquand Sun YX 


R.H. Ree, S.K. Wu 30159 (LE) 


Linzhi-Bomi, Tibet, China, alt. 3550m 


KJ808616 


KJ808645 


KJ808630 


KJ808661 


Reaumuria Linn. 


R kaschgarica Rupr. 1 


Tibet-Xinjiang Exp. Team 5166 
(HNWP) 


Tashekurgan, Xinjiang, China 


KJ808617 


KJ808646 


KJ808631 


KJ808662 


R kaschgarica Rupr. 2 


Y.M. Duan 84-A-012 (XJBI) 


Ruoqiang, Xinjiang, China, alt. 3080m 




KJ808647 


KJ808632 


KJ808663 


R soongarica (Pall.) Maxim. 


Tibet-Xinjiang Exp. Team 
5098(SHI) 


Tashekurgan, Xinjiang, China, alt. 2300m 


KJ 8086 18 


KJ808648 


KJ808633 


KJ808664 


Tamarix L. 


r. karakatensis Freyn 


K.B. Blinkovsky 12 VIII 1953 (LE) 


C. Kopetdag, Ashkhabad, Turcominia 


KJ729659 


KJ729811 


KJ 729761 


KJ729710 


r. taxo Willd. 


O.N. Demina 18 V 2001 (LE) 


Orlovsky, Bostov, Russia 


KJ 729660 




KJ729762 


KJ729711 


7. meyeri Boiss. 


M.R. Tanybaeva 12 V 2007 (LE) 


Turkestan Ridge, Kirgiztan 


KJ729661 


KJ729812 


KJ729763 


KJ729712 


7. ramosissima Ledeb. 


N.A. Brykova s.n. 10 VII 1998 
(LE) 


Orlovsky, Bostov, Russia 


KJ 729662 


KJ729813 


KJ729764 


KJ729713 



Herbaria: HNWP (Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, Qinghai, China); SHI (Shihezi University, Shihezi, Xinjiang, China), XJBI 
(Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi, Xinjiang, China); LE (Komarov Botanical Institute, Russian Academy of Sceinces, 
Russia, St. Petersburg, Russia). 
doi:l 0.1 371 /journal.pone.0097582.t001 



Estimation of phylogenetic relationships and divergence time 
was conducted using a Bayesian method implemented in BEAST 
1.5.4, employing a relaxed clock model [49-50]. We used the 
uncorrelated lognormal molecular clock model with a Yule process 
for the speciation model, GTR-I-I-I-G for the substitution model 
(estimated for the dataset), and a normal distribution with SD of 1 
as priors on the calibration nodes to accommodate for calibration 
uncertainty. Minimum ages for the two normal priors were 
constrained to the root of all taxa and a node respectively, the 
family Tamaricaceae 70 Ma, and genus Tamarix 25 Ma, with a 



detailed description as follows. A Markov chain Monte Carlo 
analysis was run for 50 mUlion generations and sampled every 
1,000 generations, and two independent runs were performed to 
confirm the convergence of the analysis. The stationarity of each 
run was examined using the effective sampling size of each 
parameter (>200). The last 40 million generations were used to 
construct the maximum clade credibility tree and the associated 
95% highest posterior density distributions around the estimated 
node ages. 
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Optimization of ancestral distributions 

Tamaricaceae root constrained. Tamaricaceae is included 
in the order Caryophyllales [2,51] and has no reliable macrofossU 
record. According to mokx ular elating [52—54], the divergence 
time of the order has been estimated as ca. 100 Ma. Tiffney [55] 
considered that the extant woody families originated during the 
Cretaceous to early Eocene, while herbaceous families appeared 
during the late Oligocene to Miocene. For instance, the woody 
families Ulmaceae and Fabaceae appeared at about 70 Ma 
[54,56] and 70-60 (-50) Ma [54,56,57] respectively. Families 
related to the Tamaricaceae, such as Polygonaceae and Amar- 
anthaceae/Chenopodiaceae, have an approximate age of ca. 
65 Ma [58]. The two subfamilies of the CaryophyUaceae have an 
approximate age of 40-55 Ma according to the age of the 
inflorescence fossil Caryophyllqflora paleogenka, and the family has an 
possible age of ca. 73 (60-80) Ma [59-60]. Even though the 
ancestor of the related families Tamaricaceae and Frankeniaceae 
has been dated to 43-30 Ma [56] , Tamaricaceae itself has had 
variable dating results. Bell et al. dated it to (72-) 60-58 (-44) [54], 
while Wikstrom et al. placed it at 52-37 Ma [56], and Schuster et 
al. at (125-) 1 18.7-1 10 (-90.7) [61]. In the light of these estimates, a 
balanced age for Tamaricaeae could be suggested as about 70 Ma; 
this estimate was chosen as the family root for molecular dating. 

The earliest reliable fossil record of Tamarix is Miocene, from 
the Yunnan province of China [62] . Most of the available fossils 
are from the Miocene, therefore, the genus might be hypothesized 
to have had an origin at least in early Miocene. However, 
considering its wide distribution in Europe, Africa, Asia, and 
North America, and our limited samples mainly from China, we 
conservatively assigned an age of late OUgocene-early Miocene, at 
ca. 25 Ma for Tamarix. 

Areas 

In accordance with the distribution of Myricaria species along the 
Asian mountains (Figure 1), we divided the distribution into six 
areas, namely. A: the eastern Himalayas, including the eastern 
QTP, the Hengduan mountains, and northern and central China; 
B: the western Himalayas, including the western QTP and the 
Pamir-Alai, Kunlun-Altun, and Hendukosh mountains; C: the 
Tianshan mountains and Junggar-Turan deserts; D: Altai-Siberia; 
E: the Mongolian Plateau; and F: Asia Minor-Caucasus-Europe. 
These six areas are distinct in biodiversity, vegetation, and 
floristics [16-18,30,32,36-38]. 

Optimization of ancestral distributions 

To infer biogeographical events, three methods were used: a 
parsimony-based procedure Diva [63], S-Diva [64] and a 
maximum likelihood-based DEC model (Lagrange; [65-66]). 
These three approaches are simultaneously considered so that to 
assess the relevant biogeographical processes, such as vicariance, 
dispersal, and extinction. 

Diva, Dispersal-vicariance analysis optimizes distributions for 
each node of the tree by minimizing the number of assumed 
dispersals and extinctions, and favors vicariance events [63,67]. 
The Diva program reconstructs widespread ancestral distributions, 
restricting them to single areas. Because allopatric speciation by 
vicariance is the null motk'l in Diva, \'icariance and range (li\'isi()n 
would always be the preferred explanation if the ancestors were 
widespread. To avoid inferring a widespread ancestor at the root 
because of the presence of widespread extant taxa, a limit of two 
areas was set (maxareas = 2) in Diva [63]. The phylogenetic 
typology' of the BEAST tree (Figure 2) was input for Diva analysis. 

S-Diva. (or Bayes-Diva) [64] is a program which comple- 
ments Diva and implements the methods of Nylander et al. [68] 



and Harris et al. [21], determining statistical support for ancestral 
range reconstructions using multiple trees from Bayesian analysis. 
This has the advantage that uncertainties in phylogenetic inference 
can be taken into account. One hundred Bayesian MCMC trees 
with the last stable typologies from BEAST, and a BEAST tree 
typology (Figure 2) were input into the S-Diva program. 

Lagrange. A valuable new biogeographical method is para- 
metric likelihood analysis, with a dispersal-extinction-cladogenesis 
model [67], as implemented in Lagrange v. 2.0.1 [65]. This 
method calculates the likelihood of biogeographical routes and 
areas occupied by the most recent common ancestor (MRCA) for 
a given phylogenetic tree topology (BEAST tree. Figure 2) and the 
present distribution of taxa. Therefore, dispersal and vicariance of 
lineages, represented by connection areas, can be estimated by the 
probabilities. This is thus a form of MRCA area reconstruction 
differing from the parsimony approach of Diva. 

Results 

Phylogenetic analysis and divergence time estimate 

The phylogenetic tree obtained from Bayesian inference in 
BEAST showed that Myricaria is monophyletic and Myrtama should 
be included in Myricaria rather than treated as a distinct genus 
(Figure 2). Within Myricaria, four clades were recognized, two 
corresponding to the existing sections Paralkhntherae Ndz. and 
Renantherae Ndz, the other two represent new groups to be named 
as sections Alpinae and Laxijlorae (see Appendbc SI). Flowers and 
filaments of the plants are illustrated in Figure 2, to show the 
characteristics of the four sections. In the present phylogenetic 
tree, the clades of the genus and the sections have strong support, 
confirming the validity of taxa at the ranks of genus and section. 
Section Renantherae comprises the most species in the genus, and 
has two subclades. The two widely distributed species of this 
section, M. alopecuroides and M. squamosa, are located in each 
subclade. The crown age of Myricaria was ca. 20 Ma, and 
8.83—6.35 Ma for four sections. 

Optimization of ancestral distributions 

The results of the three approaches Diva, S-Diva, and Lagrange 
(Figure 3) showed a consistent and strongly supported pattern, 
particularly at the ancestral nodes for Myricaria (AB), and among 
the four sections, Paralletantherae (B), Alpinae (A), Laxijlorae (A), and 
Renantherae with AB from Diva and S-Diva, whereas only with 
ABCDEF/B from Lagrange. On the whole, AB, A, and B, namely 
the Himalayas and the QTP, should be considered as ancestral 
areas in Myricaria. The events occurring in areas C, D, E, and F, 
were considered to be dispersals, several of which can be 
distinguished. M. prostrata occurs in the Himalayas, and its western 
Himalayan distribution was indicated to be a dispersal event from 
the eastern Himalayas. M. bractmta in sand lands of the Mongolian 
Plateau was shown to be a migrant from the eastern Himalayas, 
Hengduan Mountains, and Northern China. Whereas the 
distribution of M. germanica in Asia Minor-Caucasus-Europe was 
came from the western Himalayas. 

Discussion 

Phylogenetic division of sections within Myricaria 

Niedenzu [8] divided Myricaria taxa into two sections: 
Parallelantherae Ndz. and Renantherae Ndz. Gorschakova [7] 
accepted this classification system in the Flora of the USSR. 
Zhang & Zhang [4], however, considered that the establishment of 
infrageneric ranks was not appropriate due to its complicated and 
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Oligocene 



Miocene 



Pleistocene 



1 -3.77 



- Tam_ramosissima 

■ Tam_meyeri 

■ Tamjaxa 

■ Tam_karakalensis 

■ R_soongarica 

■ R_kaschgarica1 

■ R_kaschgarica2 
• M_elegans3 

- Melegansi 

■ M_elegans2 

- Mprostrata 

■ M_rosea 

■ M_wardii 
-M bracteata 



sect. Parallelantherae 



r— M_platyphylla 

*1 05 

'— Msquamosa 



J- M_alopecuroides 
'- M_paniculata 



4 -SSI 




M_alopecuroides1 
M_squamosa1 
M_germanlca 
M_pulcherrima 
M_laxiflora1 I sect. Laxiflorae 
M laxiflora2 I 





Figure 2. Chronogram of Myricaria and outgroups Tamarix and Reaumuria'm Tamariaceae, with maximum clade credibility performed 
by BEAST. Dating values are plotted at the right of the nodes, and posterior probability support of more than 95% is labeled as "*" at the nodes. 
Two vertical lines are labeled at 20 Ma and 8 Ma, corresponding respectively to two stages and two high-altitude ranges (blue ranges in A, B) of the 
QTP uplift, Himalayan motion, and rapid and major-range uplift. Four sections within Myricaria are shown, along with flowers and degree of union of 
the filaments. In section Alpinae, flower and filament status refers to M. rosea, and in section Renantherae, the flower and filaments, above right, refer 
to M. bracteata. The filaments below left refer to M. germanica and those below right refer to M. squamosa. 
doi:1 0.1 371 /journal.pone.0097582.g002 



variable morphological characters. Therefore, in the Flora of 
China [5-6] there is no division of infrageneric sections. 

However, our phylogenetic tree (Figure 2) yielded a clear 
phylogenetic division of four sections, including two that are new. 
Of them the Himalayan and QTP section Alpinae, containing Af. 
prostrata, M. rosea, and M. wardii, is characterized by the prostrate 
and recumbent, and with an adaptation to high altitudes of 3000- 
5300 m. Section Laxiflorae, comprises of only species M. taxijlora, 
endemic to the subtropical area of the Sichuan and Hubei 
provinces in eastern China. Detailed descriptions of two new 
sections are given in Appendix SI. 

Myricaria ekgans Royle was originally described from the 
Kunawar region of the western Himalayas [69]. Based on this 
species and its characters of 1 0 stamens, flat leaves, and no obvious 



style, the genus Mjrtama was established [9]. Q_aiser & Ali [70] 
named it Tamaricaria, whereas Baum [7 1] moved Myricaria elegans to 
Tamarix as T. ladachensis. As mentioned, Zhang et al. [10] and 
Gaskin et al. [1] accepted Myrtama at generic rank. Zhang et al. [10] 
considered that Myricaria elegans was an intermediate and hybrid 
genus between Myricaria and Tamarix, related more to Tamarix. 
However, Hua et al. [12] and Wang et al. [13], based on sequence 
data, found that it would be appropriately placed in Myricaria. Our 
results (Figure 1) also show that inclusion of Myricaria elegans in 
Myricaria is suitable, since the whole of Myricaria, including 
Myricaria elegans, has strong support (100%) (Figure 2). This is in 
accordance with evidence from morphological classification [4] 
and molecular phylogeny [12-14]. The former conclusion 
supporting retention of Myrtama [10] was only based on ITS 
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Figure 3. Biogeographical ancestral optimization performed with Diva, S-Diva and Lagrange. Pie charts at the internal nodes represent 
the calculated probabilities (relative frequencies) of alternative ancestral areas (reconstructions), and were produced by S-Diva; letters are labeled at 
the right and above nodes, and probabilities from S-Diva are at the right and belov*/ nodes; letters above and below at branches are from Lagrange, 
indicating the highest probability migration routes and inheritance of area by upper and lower descendant branches, respectively. Located at the 
front of M. prostrata, M. bracteata, and M. germanica, vertical lines on the branches indicate dispersals involving these three species. Area letters as 
stated in the text: A: eastern Himalayas, including the eastern QTP, Hengduan mountains, northern China, and central China of East Asia; B: western 
Himalayas, including the western QTP, and Pamir-Alai, Kunlun-Altun, and Hendukosh mountains; C: Tianshan-Tungger-Turan; D: Altai-Siberia; E: 
Mongolian Plateau; and F: Asia IVlinor-Caucasus-Europe. 
doi:1 0.1 371 /journal.pone.0097582.g003 



sequence data, which is probably not sufficient evidence [72]. 
While in the phylogeny of Gaskin et al. [1] (their Figure 2), only 
four species were sampled, and four species as clade had a strong 
bootstrap support (99%) for the inclusion in Myricaria [1]. 

Himalayan origin, ancestral inheritance, and multi- 
diversification in the Himalayas 

Our estimated crown age of ca. 20 Ma for Myricaria (Figure 2) 
falls into the probable early range of the Himalayan uplift in early 
Miocene [22], consequently allows us to speculate that uplift of this 
mountain range caused the origin of Myricaria. The biogeograph- 
ical analytical result of Diva, S-Diva and Lagrange, showing the 
combined Himalayan area AB as the ancestral area for Myricaria 
(see Figure 3), which also supports a Himalayan origin. The 
present molecular dating results are in contrast to previous 
phylogeographical opinion of a main divergence event at the 
implausible age of 1.46-2.30 Ma in the Pho-Pleistocene [14]. 

The southern and northern slopes of the Himalayas difiFer 
dramatically in temperature and precipitation. Myricaria species 
occurring on the northern slope are generally xeric, same as those 



of the main plateau [7 3] . The western portions of the Himalayas 
and adjacent QTP are more arid than the eastern parts [27,73], 
see Figure Ic. Probably these differences are the cause of the 
persistent diversification between the eastern (A) and western 
Himalayas (B) for the Myricaria lineages (Figures 1 and 3). The 
eastern A contains the alpine species M. rosea, and M. wardii, and 
the western B includes M. elegans and M. prostrata. These four 
species are endemic to the Himalayas, and occupy two of the four 
Myricaria sections. In particular, the western Himalayas (B) is an 
important geographical node and dispersal center for Myricaria, 
with movement toward the Pamir-Alai, Hendukosh, Tianshan, 
and Kunlun-Altun mountains, etc. Noticeable, the Himalayas as a 
union was divided into A and B two times (see Figure 3), once at 
the time of generic origin and diversification, and the other at the 
diversification node of the sections Renantherae and Laxiflorae. 
Overall, the Himalayan areas AB, A, and B as the ancestors 
occurred at least seven times at the nodes in Figure 3. In detail, the 
Himalayan union AB as an ancestral area appeared at two nodes, 
with estimated ages of 20.25 Ma (genus crown age) and 10.07 Ma, 
the eastern Himalayas A at three nodes with ages of 19.84 Ma, 
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6.83 Ma (section Alpinae crown age), and 6.35 Ma (section 
Laxiflorae crown age), and the western Himalayas B occupied 
two, with ages of 8.83 Ma (section Paralkkmtherae crown age) and 
5.81 Ma. All of these stressed the diversification and geographical 
heritage of Himalayan ancestry. 

Many plant groups of the Northern Temperate Zone are 
hypothesized to have originated from the QTP [16,30-36], 
especially from the Hengduan Mountains; few have originated 
from the Himalayas. Myricaria presents a case of Himalayan 
montane origin. 

In some cases, the Himalayas are regarded as a corridor of plant 
species migration linking East Asia and Central Asia and/or the 
Mediterranean [16,35,74—77]. The Himalayas have acted as a 
migration corridor for Sino-Japanese elements westward and 
Mediterranean elements eastward, such as Hippophae rhamnoides 
(Elacagnaceae) [37], Pogonophace [Astragalus, Fabaceae) [72], and 
Pkjlblohium [77] (Fabaceae). They also served as a migration route 
westwards to Central Asia from East Asia for Caragana (Fabaceae) 
[78], Begonia migration was through the Himalayas eastwards from 
Africa [40]. Our Myricaria analysis shows endemism, origin, and 
remarkable multiple diversifications within the Himalayas, but 
with weak or absent migration here. The Himalayas therefore 
seems to serve as the center of origin for Adymaria rather than as a 
corridor. This is rather like as these QTP endemic plant groups, 
such as Nannoglottis (Asteraceae) [79], Ligularia complex (Aster- 
aceae) [80], Saussurea (Asteraceae) [81], Dolomiaea, Diplazoptilon and 
Xanthopappus (Asteraceae) [82], and Aconitum (Ranunculaceae) [83], 
all of them are hypothesized to have originated from native QTP 
during Miocene-Pliocene. 

Diversification of Myricaria sections 

Our estimated crown ages for the four sections ot Myricaria were 
about 8.83~6.35 Ma (Figure 2), this is a remarkable molecular 
clock response to the rapid and major-range uplift of QTP at ca. 
8 Ma. and the diversifications of four sections all occurred in the 
Himalayas (Figure 3): sections Alpinae and Laxiflorae in the eastern 
Himalayas (A), section Paralklanthera in western Himalayas (B), and 
section Renantherae in Himalayan union (AB). We have drawn the 
distinct morphological characters of four sections in Figure 2. This 
could be regarded as the response of four sections diversifications 
to the uplift in Himalayas. 

The time of Myricaria origin and diversification related to 
Himalayan uplift in early Miocene, and the subsequent sectional 
diversifications at about 8.83~6.35 Ma, correspond respectively to 
the two labels of the QTP uplift, Himalayan motion, and rapid 
and major-range uplift [21-29]. These diversification times of the 
genus and sections are very similar to those of the genus Caragana 
[38], which also has a 16 Ma time of generic origin and 
diversification ages of sections at about 8 Ma. However, here we 
provide the Himalayan places of origin and diversification for 
Myricaria, while those of Caragana remain unknown. The evolution 
of Myricaria is also temporally similar to that of the Asian spiny 
frogs Paini [42]. Uplift of the Himalayas, is hypothesized to be 
most likely due to a cut off of the genetic exchange at ca. 1 9 Ma, 
resulting in the splitting of the subgenera jVareorana and Paa ot Paini 
began in the Miocene near 10 Ma [42]. 



Migration along Asian mountains 

For plant migration and dispersal, mountains generally act as a 
route or corridor [84] , such as the Himalayan corridor mentioned 
above. Whereas the Himalayas for Myricaria are regarded as the 
center of origin, other distributions outside of the Himalayas and 
QTP can be understood as dispersals or migrations eastward, 
westward and northward along Asian mountains (see Figure 1). In 
fact, the results of vicariance and dispersal from biogeographical 
analysis (Figure 3) show that except for divergence of phylogenet- 
ically basal clades located in the eastern and western Himalayas (A 
and B), most remaining events, occurring in areas such as the 
Tianshan-Jungger-Turan (C), Altai-Sibiria (D), Mongolian Plateau 
(E), and Asia Minor-Caucasus-Europe (F) resulted from dispersal 
events. As evidenced from Myricaria (see Figure 1), dispersal and 
migration were possible from the Himalayas to Asia Minor- 
Caucasus-Europe as shown by M. germanica, and to the sand lands 
of the Mongolian Plateau by A/, bracteata. These distributions along 
the Asian mountains are very similar to those of Hippophae 
rhamnoides (Elaeagnaceae) [37], a species with another interesting 
distribution in North Temperate Eurasia. Hippophae rhamnoides 
includes nine subspecies, and has been shown to have originated 
from the QTP, or more exactiy, the eastern QTP-Hengduan 
Mountains, and then to have radiated and dispersed in difierent 
directions. Here, Alyricaria originated in the Himalayas union, not 
in the eastern Himalayas only as H. rhamnoides. However, the 
northwestern Himalayas for H. rhamnoides and Myricaria played an 
important node role in connecting with Central Asia and Europe, 
and both dispersal route and direction is very similar. 

Supporting Information 

Appendix SI Two new sections within Myricaria. 
(DOC) 
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